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We study the tail of the spectrum for non-interacting bosons in a blue-detuned random speckle 
potential. Using an instanton approach we derive the asymptotic behavior of the density of states in 
d dimensions. The leading corrections resulting from fluctuations around the saddle point solution 
are obtained by means of the Gel’fand-Yaglom method generalized to functional determinants with 
zero modes. We find a good agreement with the results of numerical simulations in one dimension. 

The effect of weak repulsive interactions in the Lifshitz tail is also discussed. 


I. INTRODUCTION 



Effect of disorder on quantum systems has been at¬ 
tracted considerable interest in condensed matter physics 
during the last several decades. In recent years, it was 
realized that ultracold atomic gases in optical speckle po¬ 
tentials may serve as quantum simulators for diverse phe¬ 
nomena in disordered quantum systems [l|-@l- Contrary 
to the experiments in condensed matter physics, the op¬ 
tical speckles allow one to create a controllable random 
potentials acting on ultracold atoms a a. Many inter¬ 
esting features of Bose-Einstein condensates (BECs) in 
disordered speckle potentials have been addressed Jrom 
both the experimental and the theoretical side_s 
These include inhibition of transport properties 
fragmentation effects a EL El i frequency shifts ^ , 

damping of collective excitations a El El and Anderson 
localizati on Jl4H2l| . Also, the superfluid-insulator tran¬ 
sition [2214301 ] and the transport of coherent matter waves 
have been recently investigated from the theoretical point 
of view for speckles in higher dimensions [3ll436| . 

When a coherent laser light is scattered from a rough 
surface the partial waves passing through the different 
parts of the surface acquire random phase shifts. The 
interference of these randomly phased waves produces a 
speckle pattern which consists of the regions or grains of 
light intensity with random magnitude, size and position. 
The local intensity of the speckles, I{x) = |£(x)| 2 , is 
determined by the electric field £{x). To a very good 
approximation, the electric field £ (x) can be viewed as a 
complex Gaussian variable with finite correlation length 
£ giving the typical size of the light intensity grains in 
the speckle pattern. The distribution of intensity I(x ) 
across the speckle pattern follows a negative-exponential 
(or Rayleigh) law [E, [37j, 


P[I] =exp[-/// 0 ]// 0 , (1) 

where Iq = (/) is the mean intensity while the most prob¬ 
able intensity is zero. The speckle pattern shined on a 
sample of atoms creates a random potential felt by the 
atoms provided that the wavelength of the laser light is 
slightly detuned from the atomic resonance. The poten¬ 
tial is proportional to the local light intensity, V(x) = 
a J(x), so that the non-interacting atoms in the speckle 


potential can be described by the Schrodinger equation 


2 mi 


-V 2 + a I(x) 


ip(x) = E ip (x). 


( 2 ) 


where m o is the mass of atoms and ip ( x ) is the single 
particle wave function. The constant a is proportional 
to the inverse of the detuning A between the laser and 
the atomic resonance. The detuning can be either posi¬ 
tive (blue detuning) or negative (red detuning) 0. The 
blue-detuned case corresponds to a disordered potential 
consisting of a series of barriers bounded from below. The 
red-detuned speckle produces a potential bounded from 
above and made of potential wells. 

The precise single particle spectrum of the speckle po¬ 
tential is unknown even in the ID case despite the in¬ 
tense research activity in this field. It is widely believed 
that the density of states (DOS) of blue-detuned repul¬ 
sive speckles is characterized by a usual Lifshitz tail for 
potentials bounded from below [38| . However, many 
exact results known for ID random potentials [38l . l39| 
cannot be directly applied to the speckle potential since 
it is correlated and non-Gaussian. In the previous pa¬ 
per [40j ] we investigated analytically and numerically the 
ID single-particle spectrum for the both red-detuned and 
blue-detuned speckle potentials. Since the speckle pat¬ 
tern is characterized, besides the mean intensity Jo, by 
the correlation length £, this introduces a new energy 
scale E { = 7i 2 /2too£ 2 . We have shown that for dimen¬ 
sional reasons, the single-particle properties are deter¬ 
mined by the dimensionless parameter s = 2mo£ 2 a/o/J 2 
for the both red and blue detuned speckle potentials. 
We identified different Lifshitz regimes controlled by the 
dimensionless parameter s which vary from the semiclas- 
sical limit |s| 1 deep to the quantum limit |s| -C 1 [ioj . 

In the present paper we improve our results for the 
blue-detuned speckle in ID and extend them to higher 
dimensions. We apply the field-theoretic description de¬ 
veloped by Luck and Nieuwenhuizen [4lj] for studying the 
Lifshitz tails in the random potentials bounded from be¬ 
low. We generalize this approach to the blue-detuned 
speckle potential which is not only bounded from be¬ 
low but also has a finite correlation length. We argue 
that well below the corresponding energy scale Ej the 
Lifshitz tail does not depend on the precise form of the 
electric field correlations since the DOS is determined by 
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the states localized in the regions with very small inten¬ 
sity whose size is much larger than £. In this regime the 
instanton approach proposed in Ref. .41] for the poten¬ 
tial bounded from below predicts the following form of 
the DOS tail 


v(E) = A d {E) exp 


-v d F d [In (Iq/E)] {E^/Ef 2 


, (3) 


where Ud is a constant which depends on the dimension 
and F d and A d are some functions of the energy. 

The paper is organized as follows. In Sec. HJ we de¬ 
rive the replicated action for a particle in a blue-detuned 
speckle potential. The saddle point solution of the classi¬ 
cal equation in d dimensions is discussed in Sec. Mil where 
the expressions for v d and F d are calculated. The fluc¬ 
tuations around the instanton solution are investigated 
in Sec. [TV] There, the prefactor A d is calculated using 
the Gel’fand-Yaglom (GY) method [42], [43[ generalized 
by Tarlie and McKane (TM) [44| to functional determi¬ 
nants with excluded zero modes. In order to illustrate the 
power of this method, in Appendix [Cl we reconsider the 
problem of a particle in Gaussian uncorrelated disorder 
for which there are some discrepancies in the existing 
literature. In Sec. El we consider a weakly interacting 
Bose gas in a speckle potential. The obtained results are 
summarized in Conclusion. 


II. AVERAGING OVER DISORDER: 
REPLICATED ACTION 

The DOS for a particle in a particular realization of 
the electric field £{x) can be related to the imaginary 
part of the one-particle Green function 

v{E) =-Im Q(x, x; E ). (4) 

7r 

To average over different realizations of disorder potential 
we employ the replica trick ]45|. We introduce N repli¬ 
cas of the original system and use the functional integral 
representation 

Q{x,x'\E) = lim [ Vcj) <j>\{x)^>x{x') e~ s ^ (5) 
N^oJ 


with the action 

S[4>] = \ I d d xS [ ^-(^(x)) 2 + [V(x)-E}^(x)Y( 6 ) 

where 4>(x) is an TV-component scalar field. The disor¬ 
der potential is proportional to the local intensity of the 
speckles pattern created by a laser, V(x) = a £*(x)£(x), 
where £{x) is the electric field and we fix a = +1 in the 
case of a blue-detuned speckle. To a very good approxi¬ 
mation, the electric field is a random complex Gaussian 
field with zero mean and variance 


where the function G{x — y) has the width of £ and its 
precise form depends on the experimental setup [l8j, (32j . 
Then the average of the disorder potential can be ex¬ 
pressed as exponential of the sum of loop diagrams: 

exp [—y d d x £*(x)£(x)(j> 2 (x)\ 

= exp {~nO + 5iO~sO + "'}' (8) 

In Eq. ([5]) the lines stand for the correlator 0 and 
thus have two distinct ends corresponding to £* and 
£. Two lines can be connected only by ends of differ¬ 
ent types and the corresponding vertex carries a factor 
of \4> 2 {x) = \Y^n=i4 > n{ x )- The loop diagram with j 
vertices in Eq. 0 has the combinatorial factor of (j — 1)! 
which gives the number of possibilities to construct a 
loop from j lines with distinct ends. Summing up the all 
diagrams for an arbitrary function G(x) and field 4>{x) 
is a formidable task. However, there are several cases 
when one can solve this problem at least partially. For 
a special class of correlation functions G(x) the sum of 
the diagrams can be rewritten as a ratio of functional de¬ 
terminants (see Appendix 1X1). The summation of the dia¬ 
grams can be also performed within a variational method 
with Gaussian correlators and trial functions (see Ap¬ 
pendix 0 . These approximations can be used to show 
that the low-energy tail of the DOS does not depend on 
the precise form of the electric field correlator for E -C Eg 
and it is completely determined by Eg and Iq. This is in 
contrast to the Gaussian unbounded potential [45j] where 
the presence of correlations changes the low-energy Lif- 
shitz tail of the DOS (46]. For E Eg we can approxi¬ 
mate the electric field correlator by 

G(x - y) = I 0 £ d 5 d {x - y) 

where we have defined the regularized (5-function of width 
£ such that 

limS d (x - y) = 5 d (x - y) and <5|(0) = 

£—>■0 s ^ 

As a result the loop diagram with j insertions of 4> 2 (x) 
can be expressed as 


rd: 


d d x (-1)* ed 


£ d j 23 


(*>£W(z)- (9) 


i + jfWW 


By using the relation 
00 I—l V 

Ed(«W j W = -ln 

j =l d 

we can sum up all the diagrams. The averaged replicated 
action then reads 

«Sav = y d d x {^(V^(a ;)) 2 - ±E4> 2 (x) 

+-jk In [l + \(h£, d )4> 2 {x)\ | • (10) 


{£*[x)£{y)) = G(x — y), 


(7) 
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III. SADDLE POINT SOLUTION 


The field theory (flUl) has a trivial vacuum state </> = 0. 
However, a perturbative expansion around it does not 
contribute to the DOS at any finite order since the Green 
function remains real in this approximation. Following 
Ref. [4l], |45| we assume that the functional integral J5]) is 
dominated by a spherically symmetric saddle point field 
configuration. The integration over fluctuations around 
this instanton solution brings a finite imaginary part to 
the Green function Q leading to a finite DOS. 

It is convenient to express the action (1101) in terms of 
the rescaled quantities (j){x ) = \/ 2 (j)(x/£,)/^/E^ d , x = 
x/i ; and E — E/E^. By omitting the tildes on (f> and x 
we arrive at 

= j d d x | [S7 x cj>(x )] 2 - E(j) 2 (x) + hi [l + s(j) 2 (x)] | . 

( 11 ) 

The variational principle gives the following classical 
equation of motion 

V 2 0 + £^= J’- 2 . (12) 

s 1 + <fr 

Assuming that Eq. (fl2l) has a solution of the form 

4>c\(x) = n 4> 0 (x) (13) 

with n 2 = 1, we rewrite it as 

V2 * + fife = (14) 

It is instructive to compare this equation with the cor¬ 
responding saddle point equation (IC3[) in the case of 6 - 
correlated Gaussian disorder (see j45| and Appendix ICl). 
At variance with the Gaussian disorder, one cannot elim¬ 
inate the explicit energy dependence in Eq. m by any 
variable transformation. In the limit E —> 0, the classical 
solution of Eq. dl4ll approaches the form 


<A) (x) 



(s/E^/E L d (xVe) 


(15) 


in the region 0 < x V E < and essentially vanishes 
elsewhere. The functions Li(t) = cos t, L 2 (t) = J 0 (t) 
and Lz{t) = sin t/t are the spherical Bessel functions 
in d dimensions, fj, d is the first zero of L d (t*i = 7 r/ 2 , 
fi 2 = 2.40483, H 3 = 7 r) and the constants are a d = 
1, 3.71038,7T 2 in d = 1,2,3, respectively [Zlj ]. Substitut¬ 
ing this approximative solution for the saddle point into 
the action ( fTTl) leads immediately to the DOS tail of the 
form 0 , Halil 


v (E) = A d exp 


-v d (Et/E)$ ln(J 0 /£?)l, 


(16) 


with Vd = 7 r d / 2 /T(d/2 + 1). The derivation of the 

saddle point solution can be simplified in one dimension. 
The basic idea is to treat Eq. <HB as the classical action 


Scl 


dx^[V x <j) o(a;)] 2 -7/(</> 0 )} , 


(17) 


of a particle moving in the potential U((f>o) = E(j> q 2 (x) — 
In [l + s(f>o 2 (x)] with space coordinate <j >0 and time x. 
Since the system is conservative the energy of the particle 

£ 0 = 4>o 2 +U(M (18) 

is constant along any trajectory. The saddle point so¬ 
lution corresponds to the particle trajectory at £$ = 0 . 
By using a simple variable transformation the action (1171) 
can be rewritten as 

^0 

<Sci = 2 J dz\J~z~An{\A-~s~z^-~E, (19) 

0 

where corresponds to zero of the expression under the 
root. For E —► 0 the asymptotic behavior of Eq. m is 

<S c i « tt\J\/E\ n (s/Ej in agreement with the DOS (fl 6 l) 
in one dimension. 


IV. FLUCTUATIONS AROUND THE SADDLE 
POINT 


In this section we extend the instanton approach in or¬ 
der to calculate the dependence of the prefactor A d upon 
E. To this end we expand the action m around the sad¬ 
dle point (f> = (j ) c i + (j)' to second order in the fluctuation 
fields. The instanton contribution to the Green function 
is then given by 

Q{x,x'\E) ~ J Vcj)' 4> c n{x)(l) c n(x') 

xex P[-7) J d d x(/)' ol M a p(l)' p \. (20) 

Assuming the saddle point solution of the form (fl3l) the 
operator M a p can be diagonalized using the longitudinal 
and transverse projector operators in the replica space as 

A/q.^ 3 — AfL?T,Q,?T.^3 -(- AAp (S a /3 U.o.77.^) • (^1) 

The transverse and longitudinal operators can be written 
in the form 


M t ,l = -V 2 + U t ,l + 777 2 ( 22 ) 


where we have defined the mass m 2 = s — E such that 
the potentials 

s 


U T (r) = 
U L (r) = 


1 + S 4>o 
S 

1 + s (/>o 2 


— s — 


2 s 2 cj)o 


( 1+^0 2 r 


(23) 

(24) 


vanish at infinity. Note that the transverse projector op¬ 
erator has (N — 1) zero modes corresponding to invari¬ 
ance under O(N) rotations in the replica space while the 
longitudinal operator has d zero modes corresponding to 
translational invariance. In order to obtain a finite result 
from the Gaussian integration in Eq. the zero modes 
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of the operators Mt and Ml have to be separated and 
integrated out exactly without using the Gaussian ap¬ 
proximation. To that end one can perform transforma¬ 
tion to a collective coordinates Xq and n [III 0 ■ This 
yields 

G{x,x';E) ~ J d d xo dn J* (j>o{x — xo)(j)o{x' — xq) 

x J X>0'exp[-i J d d x(j)' Q M' a gij)' f3 \, (25) 

where J 4 is the Jacobian of the transformation to the col¬ 
lective coordinates xq and n and the prime in M' a p means 
that the zero modes have been omitted. The Jacobian 
calculated to leading order in the energy E by expand¬ 
ing the model in the fields around the minimum is given 
by 0 


J 1 


d d a 


(V^o) 



d d x 



(AT—1)/2 

• (26) 


The functional integral in Eq. (1251) contributes with 


V(j)'eyL p 


~2 dx< i>oM'a 


= det 'M l 1/2 • det 'M T {N 1)/2 (27) 


where det ' stands for the product of all non-zero eigenval¬ 
ues including the continuous part of the spectrum. The 
operators Mt and Ml given by Eq. (l22l) are Schrodinger- 
like operators. Their spectrum turns out to be very sen¬ 
sitive to the precise form of the saddle point solution of 
Eq. (1141) . which is not known analytically even in one 
dimension. Unfortunately, the dependence of the eigen¬ 
values on the energy parameter E cannot be easily ex¬ 
tracted using simple scaling arguments as in the case of 
Gaussian uncorrelated disorder (see Appendix ICl). Nev¬ 
ertheless, there are methods which allow one to calculate 
the functional determinants with excluded zero modes 
even without knowing precisely the spectrum. 


A. One dimensional case 

We start with the one dimensional case. The poten¬ 
tials Ut and Ul corresponding to the operators Mt and 
Ml obtained from numerical solution of the saddle point 
equation in d = 1 are shown in Fig. |Tj In the low en¬ 
ergy limit E —^ 0, we find that the potentials Ut and 
Ul approach asymptotically a square potential well of 

width ir\Jl/E and depth s =F E ^1/ln 20 respectively. 
The operator Mt has one zero mode corresponding to 
the lowest symmetric state. The operator Ml has the 
only zero energy state corresponding to the lowest anti¬ 
symmetric state while its lowest symmetric state has a 
negative eigenvalue which gives a non-zero contribution 
to the imaginary part of the Green function. 



FIG. 1. One-dimensional potential wells Ut (blue dashed) 
and Ul (red solid) for s = 0.5 and E = 0.01. 


When the spectrum of the Schroedinger operator is 
known analytically the zero mode can be explicitly ex¬ 
cluded from the product of eigenvalues. Therefore, the 
determinant can be calculated simply as an infinite prod¬ 
uct of non-zero eigenvalues. This is illustrated in Ap¬ 
pendix 1C II for the fluctuation operators Mt and Ml 
arising in the Gaussian disorder model 0 . For the blue 
detuned speckle the determinants of Mt and Ml cannot 
be calculated simply as a product of non-zero eigenvalues 
because the spectrum cannot be found analytically. For¬ 
tunately, Gel’fand and Yaglom (GY) 0 , 0 ] derived long 
ago a general formula which allows one to calculate the 
functional determinant of a Schrodinger like operator at 
least in one dimension without knowing any of its eigen¬ 
values. The GY method can be applied to an operator 
of the form 

d 2 

M =--j^+ U(x)+ m 2 , ( 28 ) 

which is defined on x € [—L, L\ for the wave functions 
satisfying the boundary conditions u(—L) = u(L) = 0. 
The limit L —> oo can be taken at the end of the calcu¬ 
lation. Since the well defined object is rather a ratio 
of two determinants than a single functional determi¬ 
nant itself it is convenient to introduce a free operator 
AJfree = — + w 2 . The GY theorem 0 ] states that 

det ^ = u{L) 

det f ree Ufree ( L ) 

where u{x) and Uf lee (x) are the solutions of the Cauchy 
problems 

Mu(x ) = 0 and Mf lee U{ re e ( x ) = 0 (30) 

with the initial conditions: 

u(-L) = Uf ree (—L) = 0, u'(-L) = Uf ree (-L) = 1. 

(31) 

Due to the presence of eigenfunctions with zero eigen¬ 
value, whose contributions to the determinant have to be 
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excluded, the GY formula (1791) has to be slightly modified 
for the operators My and M L . A simple regularization 
consists of introducing an infinitesimal shift of the spec¬ 
trum by a small shift of the mass m. Then, the original 
determinant with the excluded zero mode can be derived 
by differentiating with respect to the mass. This method 
is illustrated for the operators Mj and Ml of the Gaus¬ 
sian disorder model in Appendix 1C 21 

In the case of the speckle potential it turns out to be 
more convenient to use another regularization approach 
which has been recently proposed by Tarlie and McKane 
(TM) in Ref. [44]]. It is based on the GY method gener¬ 
alized to an arbitrary boundary conditions by Forman in 
Ref. [13] • The basic idea is to regularize the determinant 
by modifying the boundary conditions. This changes the 
zero eigenvalue to a nonzero one which can be estimated 
to lowest order in the difference between the original and 
regularized boundary conditions. Assuming that the zero 
mode of the operator A4 is given by vq(x) the ratio of the 
two determinants with excluded zero mode can be writ¬ 
ten as 


= (vq\vq) <ee(~ L ) , 3 ^ 

dctjVlfree ^o( ^free(^) 

where we defined the scalar product 

(uo|uo) := J dxvl(x). (33) 

The TM formula is very useful because the zero mode of 
the operator Mj is given by the classical solution cj>o(x) 
while for the operator Ml the zero mode is simply given 
by its derivative 4>'q{x). This is true not only in the case 
of the one dimensional speckle potential but also for the 
problem with uncorrelated Gaussian disorder where the 
TM formula is shown to reproduce the correct result of 
Cardy (see Appendix 1C 31) . By inserting the zero mode 
solutions in Eq. (1371) the ratio of the two determinants 
Mt and Ml of Eq. (1271) in one dimension can be rewritten 
as 


det' M t = {MM (j>o(-L)^{L) 

deb 1 M l l^oo {M<K >) (-L)<// 0 (L) ’ 1 ' 

because the contribution from the free operator cancels 
out. Moreover, in one dimension, the derivatives of the 
classical solution can be easily obtained from the first 
order differential equation (1151) derived from the analogy 
with the particle in a conservative potential. Although 
the function (f>o{x ) is known only numerically, the limit of 
Eq. (1371) can be calculated analytically by using that this 
solution is regular at infinity. We find the exact relation 


det' M t _ __ {MM ( _ jj} 

det ’M l {MM ^ 


(35) 


When we substitute this result in the DOS (1751) the ratio 
of two scalar products in Eq. (l35l) cancel exactly the Ja¬ 
cobian (1751) . Using the saddle point solution (fl5l) in the 



FIG. 2. The one particle DOS for a blue-detuned speckle 
potential computed numerically for s = 1 (red squares) and 
s = 0.1 (blue circles). The solid and dashed lines are the best 
fit of the tails by Eq. (1371) . 


asymptotic limit E —> 0 we obtain 


In 


Ai 


JM 

E 3/ 2 


, 1/2 


(36) 


By collecting all contributions, we find the tail of the 
DOS in one dimension 


— I 

M \EJ 


3/2 


In 


exp 


—7rs 




where A is a numerical constant and 

zo 


fi(y) = ~ J dz V 2 -1 ln(l + yz) - 1, 


(37) 


(38) 


whose asymptotics for large y is fi{y) = lny. In order 
to check this formula we have fitted the low energy tail 
of the DOS in ID speckle potential which we computed 
numerically using the exact Hamiltonian diagonalization 
in Ref. |40| . The result is shown in Fig. [7] 


B. Higher dimensions 


The GY method can be extended to determinants of 
the Schrodinger-like operators in d dimensions in the case 
of radially symmetric potentials (50j. Due to the radial 
symmetry of the operators, their eigenfunctions can be 
decomposed into a product of radial parts and hyper- 
spherical harmonics 

^ (r,#) = -irripejr) Ye($). (39) 

J' 2 

The radial parts ip#(r) are then solutions to the radial 
equations 


m t\ ^(3 ;= 


±_. 

dr 2 r 2 
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+Ut,l{t) + m 2 ) ipe(r) = A ipe(r). (40) 

The radial eigenfunctions i/j£( r ) come with a degeneracy 
factor of 


deg (£;d) 


(2£ + d-2)(£ + d-3)! 
£!(d- 2)! 


(41) 


The determinant of a radially separable operator can be 
calculated by combining the determinants for each par¬ 
tial wave £ with the weights given by the degeneracy fac¬ 
tor m as follows 


l n det^ = g deg (^) ln det, <- 


det Mf r , 


£=0 


det M, 


W 


(42) 


free 


where the free operators have been defined as 




(*+*¥) (£+*=1) 


dr 2 


m 


(43) 


In order to compute the determinants of the partial oper¬ 
ators in Eq. (143 one can use the GY method for those de¬ 
terminants that have no zero modes and the MT method 
for those that have them. The only partial operators 
which have zero modes are M and M^ 1 '. All other 
determinants can be computed using the GY formula 


det M^ l _ u w {R) 

det «&,(*) ’ 


(44) 


where and u^oe are the solution of the following 
Cauchy problems 


M w u w (r) = 0, r —>■ 0, (45) 

M free U frie( r ) = °> U free~^ + ^> r °' ( 46 ) 

The determinants of the operators M^ 0 ' and M y 1 ' with 
excluded zero modes are given by the generalization of 
Eq. (153 which reads 


iet ' M = - hm v (47) 

det Ad free R^-oo V 0 (R)U{ I ee{R) 

Here we have used that the zero modes Vo (r) of the op¬ 
erators My' and M^ 1 ' and the solution Uf lee (r ) of the 
Cauchy problem (143 have the same behavior at r — > 0. 
It turns out that the sum over £ in Eq. (143 diverges for 
d > 2 as l/£. This divergence is a g eneral property of 
functional determinants in d > 2 [5l| that was not dis¬ 
cussed in the most works on the instanton approach to 
the DOS of disordered systems Hi| HH, [52, 53j. This di¬ 
vergence reflects the fact that the field theory m has 
to be renormalized [HU that is beyond the scope of the 
present paper. To make the theory finite we introduce 
the UV cutoff A and separate the divergencies from the 
sum over £ in Eq. (143 . Following Ref. [50] we use the 
diagrammatic representation of the determinants Mt,l 
explained in Appendix O and given by Eq. (IA2I) with 



FIG. 3. The logarithm of the ratio of determinants Mt and 
Ml with subtracted diagrams A^' L computed numerically for 
s = 0.5 and d = 2 as a function of E using Eq. <138]). The solid 

blue line is a fit by 3.18156/\/I + 2.82111 In E + 10.5708. 
(Inset) The logarithm of the ratio of partial determinants for 
d — 2, E = 0.0005, and s = 0.5 as a function of £. The solid 
blue line is the asymptotic behavior 162085/(£ + 21.4092) 3 . 


V{r) = For 2 < d < 4 the only divergent di¬ 

agram is A\ which is linear in V. For d = 4 one has 
also subtract the diagram A 2 which is quadratic in V. 
Subtracting the divergent diagram from the partial de¬ 
terminants we obtain 


det Mt t , v - ^ , 

In ———— ■ — = Ai + > deg(£; d) 


det Mf r , 


£=o 


In 


det ’ M^ l 


det 


W 

free 


0(F) 

(48) 


where the symbol [-• •]0(v") means that the parts of order 
V have been subtracted. The sum over £ in Eq. (143 
becomes finite because the UV divergences have been ac¬ 
cumulated in the regularized Feynman diagram A\. The 
explicit form of the terms which have to be subtracted for 
a general potential V ( x ) is given by the expansion [5Cj | 


, detA4^ 


OO 

J drrV{r)K v (mr)I v {mr) 
0 


— J drrV{r)K 2 (mr) J dr'r' 

0 0 

xV(r')I 2 (mr') + 0(V 3 ) , (49) 

where I u and K v are the Bessel function with v = l+^ — l. 
Thus, one can compute the regularized ratio of the 
functional determinants (1481) by solving numerically the 
Cauchy problems (143- (146|) and using the GY and MT 
formulas (|44|) and (|43 . We have found that after sub¬ 
tracting the diverging term given by the first line in 
Eq. (143 the sum over £ in Eq. (143 converges asymp¬ 
totically as l/£ 3 in d = 2 and as l/£ 2 in d = 3. As an 
example, this sum is shown as a function of E for d = 2 
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and s = 0.5 in Fig. [3] The total ratio of the determinants, between the chemical potential and the density is given 

however, is dominated by the exponential of by [55| 


Af-Ai= I g [ d d x[U T {x) - U L {x)]. (50) 

J q q +m z J 

Neglecting a power-law correction resulting from the Ja¬ 
cobian (l26l) we arrive at 


A d {E) ~ exp 


Id 


h/E 
ln(/o /E) 


(rf-l)/ 2 ‘ 


(51) 


with UV cutoff-dependent coefficients 72 = 1.891n(A/s) 
and 73 = 14.81A. 


V. WEAKLY INTERACTING BOSE GAS IN 
SPECKLE POTENTIAL 


li(n) 


n 2 

2 moR 2 (n) 



(53) 


where £ = h 2 /(moC 2 ) and n c = (3 C 2 a s ) 1 . 

In the case of the speckle potential both the energy 
corresponding to the optimal well E 0 (R) = Ti 2 /{2moR 2 ) 
and the positive repulsion energy decay with growing the 
size of the typical well. Thus, there is no competition be¬ 
tween the disorder and interactions so that we have no 
need for optimization of the total energy. Neglecting the 
kinetic energy and using Eq. © to estimate the den¬ 
sity of the optimal wells we obtain with the logarithmic 
precision 


We now consider the effect of weak repulsive interac¬ 
tion on the bosons in the Lifshitz tail. We restrict our 
consideration to the three dimensional system in dilute 
regime. The corresponding Hamiltonian has the well- 
known form [55 ] 

H = J d 3 x^(x) 7^—V 2 + V(x) — fij 4’(x) 

+ |/ d d x(ipl(x)ip(x)) 2 , (52) 

where ip(x) is the secondary quantized wave function, 
V(x) - random potential and /i - the chemical poten¬ 
tial. The positive coupling constant is given by g = 
4,nh 2 a s /mo, where a is the scattering length and we as¬ 
sume a low concentration of bosons n, such that na 3 -C 1 . 
The aim is to find how the bosons fill the random poten¬ 
tial when we add the particles one by one. In particular, 
we are interested in the dependence of the chemical po¬ 
tential on the density of bosons. It is instructive to com¬ 
pare the case of a bounded from below potential with the 
random uncorrelated Gaussian potential studied by one 
of us in Ref. 155) . 

In the case of Gaussian disorder the asymptotic be¬ 
havior of the DOS for large negative energy E is domi¬ 
nated by the optimal wells of width R with the energy 
Eq(R) = —h 2 /(2moR 2 ) which decreases with shrinking 
of R. The density of optimal wells, which can be found 
from the corresponding instanton solution, is n w {R) ~ 
e~ c / R /R 3 , where C = ft 4 /(m 2 7 2 ) is the so-called Larkin 
length related to the strength of disorder 7 [55]. In the 
presence of weak repulsive interactions the positive re¬ 
pulsion energy per particle grows with decreasing R as 
E r (R) = 3gJ\f(R)/(‘iTrR 3 ), where A f(R) = n/n w (R) is 
the typical number of bosons in optimal wells. Since the 
both energies E 0 (R) and E r {R) have opposite behavior 
with respect to decreasing R one has to optimize the to¬ 
tal energy in order to find the size of the optimal well 
renormalized by interactions. This yields with the loga¬ 
rithmic precision R(n) = £/ln(n c /n). Then the relation 


g(n) « E i (V 3 In J0 7 (in ^ , (54) 

where no = ( 6 £ 2 a s ) -1 and ^3 = 47 r 4 / 3 . The asymptotic 
behavior (El holds for n no and E £ -C Io and is an 
agreement with Ref. [56j. 


VI. CONCLUSION 

We have studied the low energy behavior of the DOS 
for non-interacting bosons in a d dimensional blue de¬ 
tuned laser speckle potential. We have shown that for 
E < E% the precise form of the electric field correla¬ 
tor does not affect the asymptotic behavior. Using an 
instanton approach we have found the saddle point solu¬ 
tion which gives the leading exponential behavior. Inte¬ 
grating out the Gaussian fluctuation around this solution 
we have expressed the prefactor in the form of ratio of 
two functional determinants. In one dimension we cal¬ 
culated the ratio of functional determinants exactly us¬ 
ing the generalized GY method which allows one to take 
into account not only the discrete part of the spectrum 
of fluctuation operators but also the continuous one. In 
higher dimensions the corresponding ratio diverges that 
has been overlooked in most of the previous work on the 
instanton approach to the DOS of disordered systems. 
Using the partial wave decomposition we can separate 
the UV divergences to a regularized one-loop Feynman 
diagram and obtain a finite result for the DOS tail in 
d > 1. In Appendix 1C 41 we show that this method gives 
a correct result for the case of Gaussian uncorrelated dis¬ 
order. We also discussed the effect of weak repulsion in¬ 
teractions in the DOS tail. In contrast to the Gaussian 
unbounded disorder the interactions and disorder do not 
compete and the optimal wells are not renormalized by 
interactions that leads to a different dependence of the 
chemical potential on the bosons density. 
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Appendix A: Functional determinant representation 


The diagrams in Eq. (IA2j) can not be summed up for 
an arbitrary distribution of the electric fields including 
those that appear in real experiments. In this appendix 
we show that the sum can be performed for a special 
choice of the Gaussian distribution of the electric fields 
with zero mean and variance 

ow - / *r^u(f) 1 ’* jr ^(f)' (A1) 

where K v (x ) is the modified Bessel function. A similar 
correlator appears in the problem of the Bragg glass stud¬ 
ied in Ref. [57|]. For d = 1 the variance eb reduces to 
G(x) = Iq This is particular interesting because 

the asymptotic behavior of the DOS does not depend on 
the precise form of correlations in disorder but it is de¬ 
termined by the lower energy states which spread over 
distances larger than the disorder correlation length £. 
Therefore, in order to study the lowest order corrections 
to the asymptotic tail due to presence of correlations one 
can use Eq. (ED as a reasonable approximation for the 
variance of the electric field. 

The starting point is the following diagrammatic rep¬ 
resentation of the ratio of two functional determinants 


, ^det(-V 2 + V(x) + m 2 )\ _ (-l) n+1 A 

^ det(—V 2 + to 2 ) ) ~ ^ n n 

= CMO + iO _ - (A2) 

In the diagrams shown in Eq. E2) the dots correspond 
to the potential V ( x ) and the lines stand for the Green’s 
function 


r , x _ rn d 2 K d/2 -i{m\x\) 
0[X> (2 n ) d / 2 {m\x\) d / 2 ~ 1 ’ 

which satisfies the equation 

[-V 2 + m 2 ] C 0 (x) = S(x). 


(A3) 


(A4) 


If we separate the combinatorial factors from the dia¬ 
grams in Eq. © we obtain the same series as in Eq. (1A2I) . 
Thus, one can formally rewrite the sum of the diagrams in 
Eq. © as a ratio of two functional determinants with the 
mass to = l/£ and the potential V(x) = \Io4> 2 (x)/Cq(0 ) 
resulting from a random Gaussian electric field with zero 
mean and variance G{x) = IqCq(x) /Co(0). Note that 
divergency of Co (0) when d > 2 reflects the fact that the 


functional determinants in Eq. (IA2I) require a renormal¬ 
ization for d > 2. Here we restrict ourselves to the case 
d = 1 and obtain 


S(E) = / dx 


k (V0(a;)) 2 - ^£V> 2 (x) 


4toq 


- In 


det(—V 2 + (Jp/£) 4> 2 {x) + l/£ 2 ) 
det(-V 2 + l/£ 2 ) 


• (A5) 


In the limit £ —> 0 we can neglect the V 2 -operator in the 
determinants of Eq. (IA5I) and recover Eq. m using that 
In det = Tr In. The corresponding saddle point equation 

-^V^x) + Ej>(x) = In det [-V 2 
2?77 -q 0 (p{x) 

+ (/ o /O0 2 (*) + l/C 2 ] (A6) 


has the form of a gap equation known in relativistic quan¬ 
tum field theory |5l(. 


Appendix B: Variational method with Gaussian trial 
functions 


One can also sum up all diagrams in Eq. © for a 
particular class of functions <j>(x) that can be used to 
find variationally an approximative instanton solution by 
means of the trial function method. Let us assume that 
the one dimensional electric field correlator has the form 

G(x -y) = Io&tix - y) with d c (x) = 

For the trial function we consider 4>{x ) = h</>o(x) with 
h 2 = 1 and 

<fo {x) = VCe~™ 2/2a2 . (B2) 


By substituting the trial function into the potential part 
of the action © we find that the diagram with n electric 
field correlators reads 


^ W.2 

a 


02 = (n — 1 )! /g C n j dxi...dx n exp 
2 

7t(xi-x 2 ) 2 7r(x n _i - x n ) 2 

^ 9 ^ 

7TX 2 7r(x n - Xl) 2 


TTXi 


(B3) 


Upon making the variable rescaling Xj —> £xj/7r, the in¬ 
tegral in Eq. (IB3I) can be rewritten as 


dx\...dx n exp 


^ ) x j .4, 
i,j =1 


Xj 


jn/2 


[det A n ] 1/2 , 


(B4) 
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where we have introduced the matrix 


A 


n 


/ 2 + e 2 -1 0 

-1 2 + e 2 -1 

0 -1 2 + e 2 


0 

0 

0 


-1 

0 

0 


\ 


0 0 
V -1 0 


0 

0 


2 + e 2 -1 
-1 2 + e 2 / 


with e = £/a. The determinant of A n can then be calcu¬ 
lated, and we get 


det A n = 


m =0 


UZo(n 2 ~i 2 ) 

2 m (2 to + l)!!(m + 1)! 


2+2 m 


= 4 sinh 2 


n arcsinh- 


(B5) 


Therefore, the action evaluated using the trial func¬ 
tion (1B2I) can be written as 




S(E) = J dx 


E 


(-1)" (Io£C) r 


1 


2 n (v^) n sinh [ n arcsinh| 


• (B 6 ) 


When E ~ 1/a 2 < + = l/(+ 2 ) we can approximate 
sinh [n arcsinh^] ss ne/ 2 . As a result the second line 
of Eq. (IB 6 I) is simplified to — Li 2 [— Iq^C /This ex¬ 
pression can be also derived by applying the trial function 
method directly to the action m- This means that the 
action m obtained in the limit of uncorrelated speckle 
potential properly describes the DOS in the presence of 
correlations for E -C E^. For these low laying states the 
finite range correlations play no role because the typical 
width of the wave functions a is much larger than the 
correlation length of disorder £. 


Appendix C: The Lifshitz tail for a particle in 
Gaussian uncorrelated disorder 


In order to illustrate the power of the GY and MT 
methods, we reconsider here the problem of a particle 
in uncorrelated Gaussian disorder. There existed a dis¬ 
agreement in the literature on the preexponential factor 
in the asymptotic behavior of the DOS in the tail of the 
band 


[451 . 1531 1581 ]. The two points that have not been 
sufficiently discussed in these works are the contribution 
of the continuous part of the spectrum to the functional 
determinants and divergence of functional determinants 
in d > 2 |5§]. The replicated action of the system is given 


by |45| 


Skv = f d d x {i(V+:c )) 2 - \E(j) 2 {x) - 5 ( 7 / 4 ) 2 (</(x)) 2 } 

(Cl) 

and we look for the asymptotic behavior of the DOS in 
the limit E —oo. The saddle point solution to the 


action CD has the form 

_d = (z^E!/(V+Ba:) n, (C2) 

7 

where / satisfies the dimensionless equation of motion 
y2 /-/ = -^- (C3) 

The action (IC1I) evaluated at the saddle point behaves as 
Sd ~ (—E) 2 _d / 2 /q 2 . By repeating the calculations (l25)l - 
(1271) we find the DOS tail 


u(E) ~ (-E) 1 /2-d 2 /4+3d/4 


I det 7 Mt 
det 7 Ml 


e~ Scl , (C4) 


where the regularized determinants of the fluctuation op¬ 
erators need to be calculated. The transverse and longi¬ 
tudinal operators derived by expansion around the saddle 
point solution have the form 



(C5) 


(C6) 


where the energy E is assumed to be large and negative. 


1. Brute force method in d = 1 

The one dimensional case is interesting for testing the 
GY method because the ratio of determinants of the op¬ 
erators (1C5I) and (1C6I) can also be calculated directly from 
the product of their eigenvalues. The saddle point solu¬ 
tion of Eq. (IC3I) is 


f(x) = 4-\/3 sech(a;) (C7) 


and the operators (IC5I) - (IC6I) can be rewritten in the form 
of Poschl-Teller operators [43l. [HI 




— j (j + 1 ) sech 2 (^/—E x'j 


(C8) 


where j takes integer values. The operators (IC5I) and 
(|C 6 l) correspond to the case Mt = Alia and Ml = 
A 4 i, 2 - The spectrum of the Poschl-Teller operator A i m j 
contains a discrete part which have j bound states 
£( = \E | (m 2 — £ 2 ) (£ = 1 ,..., j) and the continuous part 
e(k) = \E | (fc 2 + m 2 ) with the density of states 


i k ) 


L 

2ir 


1 

7r 


E 


t=i 


i 

P + k 2 ’ 


(C9) 
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which has been regularized by putting the system in a 
box of size L. This yields 


In order to exclude the zero modes, we apply a shift of 
the mass m —> Vm 2 + Sm 2 that gives 




1=1 


/ oo 

dk Vj (k) log [|.E| (fc 2 + m 2 )], (CIO) 

-oo 


where the free particle contribution L/2tt has been can¬ 
celed. From this formula the functional determinants of 
the operators (|C5l) - (|C6l) can be determined straightfor¬ 
wardly. The operator Mt = Mis lias just one discrete 
zero eigenvalue. Excluding this zero mode, the only con¬ 
tribution comes from the continuous spectrum given by 
the second line in Eq. (IC10I) . 


det {Mm,j + Sm 2 ) 2 det 'M m ,j 
det {M m fi + Sm 2 ) det M m ,o ’ 

Restoring the dependence on E we obtain 


det'.Mi j 
detM^o 


(- 1) J ' +1 


1 

2 j{j + i)\e\ 


(C17) 


in agreement with Eq. (IC11I) - (IC12I) for j = 1 and j = 2 
respectively. 


3. McKane-Tarlie formula 


det 'Mt 
det Mf ree 




(Cll) 


where Mf ree := M m , o- The operator Ml = M\^ has two 
discrete eigenvalues: a negative eigenvalue (£ = 2) giving 
a finite imaginary part of the Green function (U), and 
a zero eigenvalue (£ = 1). The negative mode and the 
continuous spectrum contribute with 

1 


det 'M l __ 

det Mfree 12|i£|' 


We obtain 


lim det 'Ml 1/2 • det 'M T {N 1)/2 

iV->0 

which is independent of E. 


2. Regularized Gel’fand-Yaglom formula 


We now apply the MT method. First, we need the 
zero modes of the operators Mt = Mp and Ml = 
M\s- They are given by v\s{x) = \E\f(V~Ex) and 
^ 1 , 2 ( 2 ;) = \E\^ 2 f'{V~Ex), respectively. Then, accord¬ 
ing to Eq. (l34l) . the ratio of the determinants with ex¬ 
cluded zero modes is given by 

det'M T _ (vi,i|ui,i) v' 12 {-L)v' 12 {L) o 
clot'M L " <wi, 2 |wi, 2 > v[ i(-LK 1 (L)~ ’ 


(C12) 

where we used 

(CIS) 

i Vs, (C13) 

(wi,i|vi,i) = 2 |E| 3/2 , 

O 

(C19) 


{Vl,2\Vl, 2 ) = g|-B| 5/2 , 

(C20) 


^1,2 ( OO) _ V' 12 (00) 

v l,l(-°°) <,i(oo) 

(C21) 


The same result can also be obtained using the GY 
method (l29l) . The solution u m j of the corresponding 
Cauchy problem m-m for the Poschl-Teller opera¬ 
tors M m ,j (1C8I) can be found analytically. For the sake 
of compactness, we show here the formula for |i?| = 1, 

U m ,j{L) = tt : : 7 

{l+j-m) 

-If (Y) Qf (-Y) + P™ (-Y) Qf (Y) 

X P™ +1 i-Y) Qf (-Y) - If (—Y) Qf +1 {-Y) ’ 

(C14) 


where we defined Y := tanhL and Ff{x) and Qf(x) are 
the Legendre functions. The solution (|C14j) for the free 
operator A4f ree = M m ,0 reduces to 


sinh [2mL\ 
u m ,o{L) = - 

m 


(C15) 


The ratio of the two determinants is then given by 


det ~+i Tn j _ Um,j{L) 

detA4 mj o l-j-oo u m 0 (L) 


(C16) 


4. Gel’fand-Yaglom method generalized to radial 
operators for d > 1 

The radial parts of the eigenfunctions of the opera¬ 
tors (IC5I) and (IC6I) satisfy Eq. (1401) with the mass m=\E\ 
and the potentials 

U T {r) = (/=Er)\ (C22) 

U L {r) = ^f(V^Er ) 2 . (C23) 

Scaling analysis shows that the solutions of the corre¬ 
sponding Cauchy problems (1451) and (1461) have the form 

u w = \E\g^{\E\^ 2 r), (C24) 

^L = \E\9^ e (\ E \ 1/2 r)- (C25) 

There is no zero modes for £ > 1 and we can apply the 
GY formula (l44l) . The ratio of the partial determinants 
for £ > 1 is given by 

lim gW(\E\ 1 ' 2 R)/gg e (\E\ 1 / 2 R), (C26) 

R—y 00 
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which does not depend on E. Thus, all the ratios of 
the partial determinants with £ > 1 do not contribute 
to the E dependance of the full ratio of the functional 
determinants. 

The operators M^ and have a zero eigenvalue 
so that to exclude it we apply the MT method. The 
corresponding ratios of the partial determinants can be 
calculated using Eq. S3)- The scaling behavior of the 
zero mode eigenfunction vo ( r ) is again given by 

v o( r ) = |£%o(|£| 1/2 r). (C27) 


by Eg. (Id9l) and read 

OO 

j dr rU T Ar)K v {\E\ 1/2 r)I,{\E\ i,2 r) (C32) 

o 

where v = l + | — 1 and Ut,l are given by Eqs. (1C22I) 
and (IC23I) . It is easy to see that the terms (IC32I) do not 
depend on E. The bare diagram Ai can be written as 

Af’ L = [ d d xU T ,L{\x\) lim C 0 (x - y), (C33) 

J y^ x 


This yields 

(v 0 \v 0 ) ~ \E\ 3 / 2 , (C28) 

lim v' 0 (R)u hee {R) = \E\ 5 / 2 lim g' 0 (R)g hee {R),( C29) 

R—too R —^oo 


where the last limit is expected to be finite. Thus, the 
logarithms of the determinant ratios are given (up to an 
energy independent constants) by 


1 det 'Mt 

det Ibf free 


deg (0; d) 


, det'M® 


In |J5| _1 

(C30) 


and 


det 'M l 
In -——— 
det A/f ree 


deg (1; d) 


, det ' 

ln - hr 


dlnlE^ 1 . 

(C31) 


Above we assumed that the ratios of the determinants 
are finite and the sum over £ is converging. However, we 
know that this sum diverges for d = 2 and d = 3. We have 
to subtract from each ratio of the partial determinants 
the term resulting from the partial wave decomposition 
of the diverging diagram Ai. After that the regularized 
diagram Ai has to be added to the action as shown in 
Eq. (H51) . The terms needed to be subtracted are given 


where Cq(x) is given by Eq. (IA3I) . The expression (1C33I) 
diverges and has to be regularized, e.g. by the UV cutoff 
A as follows 

A? L =[ - 2 TZT 2 f d d xU T , L (\x\), (C34) 

J\q\<A H + m J 

where the last integral behaves as (— E) l ~ d / 2 . Combining 
Eqs. (IC30I) - (1C31I) and (1C34I) we obtain 


det 'Mt 
det 'M l 


i\E\ — e 


(Af-Af )/2 


(C35) 


where the factor i comes from the negative eigenvalue 
of the partial operator m|°\ Inserting Eq. (1C35I) into 
Eq. (IC4I) gives the DOS tail for E —> —oo. Additionally 
to the term {—E) 2 ~ d ^ 2 resulting from the action evalu¬ 
ated at the saddle point the exponential now contains the 
term (— E) 1 ~ d ^ 2 with a non-universal coefficient which 
depends on the UV cutoff. Fortunately, we already know 
how to renormalize the field theory CD which is nothing 
but </> 4 theory. It easy to recognize in the first integral 
in Eq. (|C34|I the one-loop diagram which shifts the mass 
to 2 = \E\. The diagram (IC34I) is compensated by the 
counterterm coming from the mass shift. Thus, in terms 
of the renormalized energy Er = E — Eq, where Eq is 
some non-universal energy scale depending of the UV 
cutoff, the DOS tail has the form 


v(Er) ~ (-E R ) d< - 5—0/4 e - c onst(-B K ) 2 d/2 /7 2 . ( C3 g) 
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